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PURPOSE 


The  purpose  of  the  project  is  to  increase  understanding  of  the  roles  of 
local  circulations  in  shaping  patterns  of  atmospheric  fallout  at  the  earth’s 
surface. 


ABSTRACT 


Model  local  circulations  are  discussed  with  regard  to  their  influence  on 
the  disposition  of  airborne  particles  in  a  non-precipitating,  non-diffitsive 
atmosphere.  Two-dimensional  vertical  circulations  like  the  sea  breeze  along 
a  straight  coast,  or  like  the  circulations  near  mountain  ranges,  when  steady- 
state,  principally  affect  the  time  of  establishment  of  a  pattern  and  not  its  shape 
or  the  ultimate  patterns  of  the  total  accumulations  at  the  ground.  When  such 
calculations  intensify  or  decay  during  the  passage  of  fallout  through  them,  both 
the  duration  and  onset  time  of  fallout  at  ground  locations  can  be  affected,  with 
associated  local  changes  of  total  accumulation. 

Three-dimensional  local  circulations  affect  the  place  of  deposit  of  air¬ 
borne  particles  and  the  time  of  their  deposition.  Circulations  confined  to  the 
horizontal  plane  can  affect  the  place  of  deposit  and  shape  of  a  fallout  pattern, 
but  vertical  air  currents  must  occur  to  change  the  time  of  its  establishment. 

Three-dimensional  divergence  accompanying  vertical  currents  in  the  real 
compressible  atmosphere  significantly  increases  local  concentrations  in  only 
those  cases  where  updrafts  are  about  equal  to  particle  fall  speeds;  in  such  cases, 
the  role  of  diffusion  should  be  simultaneously  considered. 

The  influence  of  local  circulations,  in  any  case,  is  greatest  with  respect 
to  the  smallest  (slowest  falling)  particles.  C  ) 
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1.0  FACTUAL  DATA 


1.1  Introduction 

The  distributions  of  close-in  fallout  resulting  from  atomic  explosions  have 
usually  been  considered  from  the  viewpoint  of  the  descent  of  size-distributed 
particles  in  the  temporally  uniform  wind  field  defined  by  radiowinS  reports 
[1, 2,  3, 4,  5].  Effects  of  widespread  and  localized  precipitation  have  been  con¬ 
sidered  by  some  investigators,  notably  Austin  [6]  and  Greenfield  [7].  Although 
the  importance  of  local  wind  circulations  has  been  recognized,  especially  in  the 
case  of  high-yield  bursts,  there  has  been  little  treatment  of  kinematic  processes 
fundamental  to  the  effects  of  local  circulations  on  fallout  distributions.  The 
present  study  is  intended  to  fill  this  gap,  and  analyzes  model  circulations  to  pro¬ 
vide  understanding  essential  for  prediction  of  the  effects  of  local  circulations 
and  fallout  in  the  real  world. 

Like  Kellogg  et  al.  [3] ,  consider  the  distribution  of  contaminants  in  a 
stabilized  cloud  that  exists  as  stem  and  mushroom  some  minutes  after  a  nuclear 
detonation.  The  contamination  products  are  assumed  to  move  horizontally  with 
the  environmental  wind  u,  while  failing  at  speed  relative  to  the  air.  V^,  a 
function  of  height,  is  a  negative  number  when  it  refers  to  the  fall  velocity. 

Figure  1  represents  a  portion  of  the  stem  cloud  within  which  particles  of 
various  diameters  and  fall  speeds  are  distributed.  If,  for  example,  the  wind 
speed  increases  at  a  constant  rate  with  height,  the  trajectories  of  particles  of 
fall  speed  V  starting  at  the  three  points  A^,  B^,  are  like  the  parabolic  arcs 
shown.  Since  the  environmental  wind  field  is  assumed  steady,  all  particles  of 
fall  speed  V.  following  trajectories  which  pass  through  points  A  ,  B  ,  and  C  , 

1  6  w  6 

respectively,  keep  to  the  same  respective  paths  beyond  the  location  of  the  initial 
cloud.  It  is  appropriate  then  to  consider  a  suitably  large  number  of  points  on 
the  leeward  side  of  the  stem  cloud  at  the  initial  time;  and,  for  each  point,  the 
corresponding  mass  6M.  of  fallout  material  which  falls  at  speed  and  is  con¬ 
tained  in  a  stream  tube  of  unit  horizontal  cross  section,  extending  from  the  lee¬ 
ward  points  to  the  windward  boundary  or  top  of  the  cloud.  Since  the  wind  field 
in  this  first  instance  is  nondivergent,  the  stream  tube  of  unit  horizontal-cross¬ 
section  passes  through  each  point  with  unaltered  cross  section  to  the  ground. 
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Idealized 
shape  of 
stem  cloud 


The  mass  per  unit  area  associated  with  each  point  can  be  represented  by  the 


equation 


6M.  *  fC.  sin  ©  ds  =  fC.  dz 


where  C.  is  the  concentration  of  material  of  fall  speed  V.  at  each  point  along  the 

1  1  2  2  1/2 
trajectory  within  the  cloud  at  the  initial  time,  sin  0  =  -V./ (V.  +  u  )  ,  ds  is 

incremental  arc  length  along  the  trajectory,  and  dz  is  the  height  increment. 

Subscript  s  denotes  that  the  integration  takes  place  along  the  trajectory. 

For  each  point,  there  is  a  corresponding  mass  of  material,  and  we  now  have 

to  define  when  and  where  it  reaches  the  ground.  Obviously,  the  first  particles 

along  a  particular  trajectory  in  the  figure  reach  the  ground  when 

first  j0  V’  {) 


and  the  last  when 


f  1  dz 

0  v{ 


In  Fig.  1,  where  Z  and  Z  are  the  heights  at  which  the  trajectory  inter- 

£*  A 

sects  respectively  the  right  and  left  hand  boundaries  of  the  atomic  cloud.  The 
duration  of  fallout  is 

z 

_  f  1  dz 

T  =  -  /  —  | 

duration  zn  V. 

0  i 


The  negative  signs  in  Eqs.  (2),  (3),  and  (4)  correspond  to  our  practice  of 
representing  the  fall  velocity  V  by  a  negative  number.  While  the  particles  are 
descending,  they  are  displaced  to  a  horizontal  distance  X  by  the  horizontal 


wind: 


X  =  /T  udt  =  f  -7-  • 


z  V. 
1 


Therefore,  when  the  horizontal  winds  within  the  stabilized  atomic  cloud  are 
defined,  any  point  on  the  cloud  edge  can  be  associated  with  a  mass  per  unit 
area  of  fallout  of  fall  speeds  near  V  ;  a  point  on  the  ground;  and  (since  the 
heights  of  the  end  points  of  the  trajectory  within  the  atomic  cloud  are  defined) 
the  times  of  termination  of  fallout  at  the  ground  (see  Fig.  1).  Isopleths  drawn 
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to  fit  the  values  attached  to  points  at  the  ground  define  the  pattern  of  fallout  due 
to  particles  falling  in  the  range  of  speeds  about  V^.  The  addition  of  patterns 
pertaining  to  calculations  for  particles  of  different  fall  speeds  furnishes  the  time- 
dependent  pattern  of  total  fallout  concentration  at  the  ground.  The  isopleths  of 
mass  can  be  converted  to  levels  of  radioactivity  from  a  relationship  expressing 
the  radioactivity  of  unit  mass,  of  the  form 

R  =  RQe~kt ,  (6) 

where  R^  and  k  depend  on  the  nature  of  the  nuclear  detonation. 

Methods  developed  by  various  investigators  differ  somewhat  in  their 
detailed  assumptions  and  lead  to  forecast  fallout  patterns  which  are  different 
from  each  other  and  from  observations  in  the  real  world  14 J.  Consideration  of 
the  effects  of  local  circulations  may  suggest  means  for  improving  the  existing 
fallout  models. 

This  report  discusses  how  local  circulations  outside  the  atomic  cloud  can 
alter  the  mass  of  fallout,  the  time  of  onset,  and  points  on  the  ground  associated 
with  points  on  the  boundary  of  the  stabilized  cloud.  The  treatment  assumes  as 
above  that  the  diffusion  is  unimportant,  though  cases  occur  in  which  the  omission 
of  diffusion  must  lead  to  relatively  large  errors.  The  study  uses  idealized 
models  to  facilitate  mathematical  treatment  and  the  development  of  physical 
understanding. 

1.2  Continuity  Equations  for  Airborne  Particles 

Consider  air  moving  at  a  velocity  having  components  in  the  x,  y,  and  z 
directions  of  u,  v,  and  w  respectively  (w  is  positive  in  the  upward  direction).  The 
air  contains  particles  of  concentration  C,  which  fall  relative  to  the  air  at  the 
velocity  ,  a  negative  quantity.  Then  a  development  similar  to  Haurwitz  [8] 
gives  the  following  fundamental  continuity  equation  for  the  concentration  of 
particles  during  their  descent  towards  the  surface. 

<c<» +  v,)»  ™ 

Equation  (7a)  assumes  that  particles  whose  concentration  is  C  are  neither 
created  nor  destroyed  on  their  way  down  and  that  they  are  distributed  homogeneously 
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over  a  volume  of  air.  A  similar  continuity  equation  holds  for  the  motion  of  the 
air 

a  =  -  i  <*■>  *  ^  W  *  &  «»•>•  <7b> 

In  the  atmosphere  it  is  practically  always  true  that  horizontal  changes  and  local 
time  changes  of  the  air  density  are  relatively  small.  By  omitting  these  terms, 

Eq.  (7b)  can  be  simplified  to 


cfci  dv  dw  3  In  p 

—  +  —  +  —  =  -  w  -  . 

dx  dy  dz  3z 

Substitution  of  (7b)  into  (7a)  yields 


c.f.  Kessler  [9].  The  derivatives  of  the  wind  appearing  in  Eq.  (7a)  are  absent 
from  Eq.  (9);  the  density  term  which  replaces  them  accounts  for  the  compressi¬ 
bility  of  the  atmosphere. 

Equation  (9)  defines  the  change  in  concentration  at  a  fixed  place  in  the  air 
in  terms  of  the  winds  and  fall  speeds  at  the  place  and  the  gradients  of  the  con¬ 
centration  near  the  place.  For  many  purposes,  however,  it  is  convenient  to  be 
able  to  refer  to  changes  of  concentration  following  the  motion  of  particles.  A 
local  change  can  always  be  equated  to  the  sum  of  individual  and  advective  changes, 
i.e., 


ac  dc  ac  ac  ,  „  v  ac 

at  dt  3x  3y  '  r  9z 
Substitution  of  Eq.  (10)  into  Eq.  (9)  yields 

av 

dlnC  i  alnp 

-  _  _  -  +  w  - . 

dt  az  az 


(10) 


a  i) 


This  equation  is  a  cornerstone.  It  says  that  only  two  factors  (diffusion 
omitted)  affect  the  concentration  of  the  set  of  fallout  particles  falling  at  speed 
V..  One  factor  is  the  vertical  variation  of  the  fall  velocity  V.;  the  other  is  the 

l  l 

atmosphere’s  compressibility  (i.e.,  three-dimensional  divergence  of  vertically 
moving  air).  Since  the  fall  velocity  of  a  particle  varies  principally  as  a  result 
of  its  entering  denser  air,  both  terms  in  Eq.  (11)  can  be  related  to  the  compres¬ 
sibility. 
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1.3  The  Effects  of  Steady  Vertical  Air  Circulations  in  an  Incompressible 
Atmosphere 

Since  the  terms  on  the  right  side  of  Eq.  (11)  usually  have  magnitudes  of 
-4  -1 

the  order  10  sec  ,  the  concentration  of  airborne  fallout  is  not  importantly 
affected  by  them  over  the  course  of  an  hour  or  so.  In  an  incompressible  atmos¬ 
phere,  the  fall  speed  of  a  particle  relative  to  the  air  is  constant  and  Eq.  (11) 
becomes 


dlnC 

dt 


=  0; 


C  =  constant. 


(12) 


Equation  (12)  is  accurate  for  circulations  whose  depths  are  a  few  kilometers  or 
less,  and  where  vertical  currents  are  considerably  smaller  than  the  fall  speeds 
of  the  particles. 


1,3.1  Descent  of  Particles  Through  Vertical  Air  Currents  Without 
Horizontal  Motion 

It  is  obvious  that  particles  falling  through  steady  vertical  currents,  where 
the  horizontal  air  motion  is  zero,  arrive  at  the  ground  at  the  same  location  that 
would  have  been  their  destination  in  the  absence  of  air  currents.  The  effect  of 
vertical  currents  on  fallout  patterns  is  to  change  the  time  when  they  are  estab¬ 
lished.  The  time  T^  for  particles  with  fall  velocity  V.  to  descend  through  a 
layer  of  depth  H  without  a  vertical  circulation  is  T^  =  -(H/V.).  When  there  is 
a  vertical  circulation,  the  time  is  approximately  T  =  -[H/(V.  +  W)J. 

If  the  updraft  w  in  the  vertical  column  anywhere  exceeds  V.,  particles  can¬ 
not  reach  the  ground.  When  updrafts  exceed  fall  speeds,  however,  other  factors 
become  important.  These  include  diffusion,  changes  of  the  horizontal  wind  field 
which  carry  the  particles  out  of  the  updraft  column,  and  the  effects  of  three- 
dimensional  divergence  due  to  the  compressibility  of  the  atmosphere. 

The  time  taken  for  particles  of  various  fall  speeds  to  fall  through  vertical 
currents  defined  by  (see  Kessler  {9 J) 


dz 

dt 


4w  2 

max  .  z  . 

=  *  -  -jj-  <l  -  TT> 


(13) 


can  be  found  by  integrating  the  sum  of  this  expression  with  the  fall  speed  between 
the  upper  boundaries  of  the  circulation  and  the  surface 
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t  m 


0  |4w  2 

,  I  max  z 

'H  pr-(I  -  \ 


d z.  (14) 

The  time  taken  for  a  particle  to  travel  through  the  updraft  core  of  the  circulation 


is 


H 


up 


1/2 


w  (-  — 
max  w 


-  1) 


tan 


-1 


<- 


max 


w 


l/2j 

-  1)  ^ 


(15) 


max 


The  travel  time  through  the  downdraft  is  given  by 


H 


1  + 


(1  - 


down 


1/2 


In 


jCI 

|w|  '  J 

max 


|w| 


max 


4|w|  (- 

max  |w| 


♦  1) 


max 


1  - 

. 


(1  - 


v.  1/5r 

— i — j 
|wi 

max  J 


w 


max 


in  which  |w|  is  the  maximum  absolute  value  of  the  vertical  velocity, 
max 

Figure  2  shows  the  travel  times  in  relation  to  the  particle  fall  speeds  and 
the  vertical  velocities  of  the  air. 


(16) 


1.3.2  Role  of  Complete  Cellular  Circulations  with  Vertical  Air  Currents 
Of  course  the  presence  of  vertical  currents  implies  the  existence  of  hori¬ 
zontal  divergence  and  horizontal  air  velocities.  The  updraft  ana  downdraft  columns 
treated  in  the  previous  section  occur  at  particular  horizontal  positions  in  a  circu¬ 
lation  depicted  in  Fig.  3.  If  defines  the  maximum  vertical  velocity,  z  =  0 

and  z  =  H  the  lower  and  upper  boundaries  of  this  two-dimensional  circulation 
and  L  its  width,  the  conditions: 

a)  da/d\  +  9w/9z  =  0 

b)  w  =  0  for  z  =  0  and  z  =  H 

c)  w  =  w  for  z  =  H/2 

max 

are  satisfied  by  equations  of  the  following  type 


u 


w 


2Lw 

max 

iH 

4w 


1)  sin 


2 


max 

H 


2  TX 
L 

2  TX 
L 


(17a) 

(18a) 
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for  the  motion  of  air. 

With  addition  of  an  ambient  windfield  =  f(z),  Eq.  (17a)  becomes 
.  2Lw 

dx  max  ,2z  ...  2  irx  .  ...... 

ST  =  “liT-  (_S  ' 1J  9U1  ~T  *  |f(z)|-  (17b) 

For  particles  descending  through  the  circulation  with  a  fall  speed  V  ,  Eq.  (18a) 
becomes 

4w  2 

dz  max  z  .  2  irx  . 

*  -  ~T~ <z  -  H> cos  *  vr  <18b) 

It  is  possible  to  represent  the  trajectories  of  the  particles  by  a  stream  function 
ip,  which  satisfies  8u/8x  +  8(w  +  V.)/3z  =  0,  by  putting 


jto  dip 

at  =  '  az 


and  (w  +  V.) 

l 


M 


If  this  stream  function  indeed  exists,  the  following  equation  must  be  satisfied: 


^  =  -  /  —  dz  =  /(w  +  V. )  dx. 


(19a) 


_  2Lw 

,  8x  ,  max 

-  f  —  dz  =  - - 

J  at  ttH 


M  - 

(z  -  ~)  sin  ~~  -  /f(z)  dz  +  f(z) 


(19b) 


/(w  +  V.)dz 


4w  2 

max  .  z  .  L  2  nx 

IS-  (z  -  TT>  T*  sm  IT  *  V  +  t<z)- 


(19c) 


Equations  (19b)  and  (19c)  are  equal  if 


2Lw  2 

max  .  z  .  2  irx  rc,  ,  , 

^  -  (z  -  — )  sin  —  +  V.x  -  Jf(z)  dz.  (19d) 

In  this  case,  the  fall  velocity  V  is  assumed  to  be  constant.  However,  if  V.  is  a 

function  of  z  it  will  not  be  possible  to  derive  ip. 

It  is  instructive  to  plot  the  streamlines  of  fallout  descending  at  constant 

speed  Vj  through  the  circulation  defined  by  Eqs.  (17b),  (18b),  and  (19d)  with 

f(z)  =  0.  Figure  4  shows  the  streamlines  when  -(V./w  )  =  2.  Note  that  each 

i  max 

streamline  of  fallout  terminates  at  a  point  directly  beneath  the  point  where  it 
enters  the  circulation.  In  other  words,  this  local  circulation  does  not  affect  the 
pattern  of  concentration  at  the  ground,  but  only  the  time  of  its  establishment. 
Sediment  descending  through  updrafts  arrives  later  than  would  otherwise  be  the 
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Horizontal  distance  x,  km 

Fig.  4.  Trajectories  of  fallout  descending  through  the  local  circulation  pictured  in  Fig.  3,  where  -V  /w 


case,  and  that  descending  through  downdrafts  arrives  earlier.  At  x  =  I./2  and 
x  =  L  (see  upper  boundary  of  Fig.  4),  the  arrival  times  are  given  by  Fig.  2, 
and  elsewhere  there  is  a  smooth  variation  between  these  extremes. 

That  the  locations  of  the  exit  points  of  fallout  streamlines  are  unchanged 
by  any  steady  local  circulations  where  V.  =  0  and  the  incompressible  approxima¬ 
tion  applies  (i.e.,  cases  described  by  Eqs.  (17b),  (18b),  and  (19d)]  is  readily  seen 
by  considering  implications  of  the  converse.  If  the  exit  points  were  changed,  the 
termini  of  some  of  the  streamlines  would  be  comparatively  crowded;  elsewhere 
comparatively  separated.  This  implies  a  change  of  the  concentration  C.  But 
Eq.  (12)  says  that  C  remains  unchanged.  Figures  5a,  5b,  6a,  and  6b  show  that 
the  entrances  and  exit  points  of  streamlines  of  fallout  are  also  unchanged  by  the 
presence  of  a  local  circulation  in  an  atmosphere  with  ambient  shear. 

The  invariance  of  the  fallout  pattern  can  also  be  proved  mathematically 
as  follows.  The  local  wind  field  under  consideration  satisfies  the  continuity 
condition 
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Fig.  5a.  Schematic  representation  of  streamlines  in  a  two-dimensional 
shearing  flow'  without  a  local  circulation. 
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Fig.  5b.  Trajectories  of  fallout  descending  through  the  wind  field 
depicted  in  Fig.  5a. 
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Fig.  6a.  Streamlines  simulating  a  kind  of  airflow  over  a  mountain  in  an 
atmosphere  with  wind  shear. 
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Fig.  6b.  Trajectories  of  fallout  descending  through  the  circulation 
depicted  in  Fig.  6a. 
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Also,  >p(x,z)  must  be  constant  on  the  boundary  of  R  because  of  conditions 
Eq.  (21),  i.e.,  if  (x^z^)  and  (x^z^)  are  both  on  boundary  of  R,  ^(x^z^)  * 
fl^z  ).  It  follows  that  the  general  integral  of  Eq.  (22)  is 

zi 

ip(x,  z)  +  V.x  -  /  f(s)ds  .=  C  .  (24) 

1  z0  1 

Hence  a  particle  that  enters  R  at  the  point  (xqzq)  describes  the  curve 
z 

V.x  -  /  f(s)ds  =  V.x  .  (25) 

■i  zQ  »  0 

Equation  (25)  says  that  the  exit  point  depends  only  on  the  entrance  point,  the 

fall  speed  V.,  and  the  environmental  wind  field  f(s).  In  particular,  it  is  independent 

of  the  local  wind  field. 

The  important  implication  of  the  above  is  that  such  idealized  steady-state 
circulations  as  the  two-dimensional  mountain  wave  model  and  the  sea  breeze 
parallel  to  a  straight  coast  have,  in  the  absence  of  diffusion,  only  the  effect  of 
altering  the  time  of  establishment  of  a  fallout  pattern,  and  have  no  effect  on  the 
ultimate  shape  of  the  pattern  or  its  concentration.  It  is  shown  in  Section  1.5 
that  the  variable  fall  velocity  of  particles  and  the  compressible  nature  of  the 
real  atmosphere  only  slightly  alter  this  conclusion. 

Such  a  simple  theorem  does  not  hold  for  a  steady  axially-symmetric 
vertical  circulation;  e.g.,  one  with  updrafts  in  a  central  region  and  a  surrounding 
ring  of  subsiding  air  which  is  embedded  in  a  general  flow.  An  example  of  such 
a  circulation  is  shown  in  Fig.  8a.  Figure  8b  illustrates  the  flow  resulting  from 
the  addition  of  the  circulation  in  Fig.  8a  to  a  uniform  flow.  A  particle  at  point 
P  in  the  x  -  z  plane  (Fig.  8b)  follows  the  trajectory  PQT  and  will  remain  in  the 
x  -  z  plane  during  its  descent.  A  particle  arriving  at  R,  however,  will  not 
remain  in  the  plane  ABCD,  (parallel  to  the  x,  z  plane)  since  the  local  circulation 
carries  it  to  a  point  S  wherefrom  it  follows  the  direction  of  the  ambient  flow  in 
the  x-direction.  It  therefore  does  not  arrive  at  that  point  on  the  ground  that 
would  be  its  destination  if  the  local  circulation  were  zero.  This  illustrated 
circulation  thus  redistributes  the  particles.  When  such  circulations  are  large 
compared  to  the  extent  of  fallout,  the  pattern  of  fallout  may  be  significantly 
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Fig.  8a.  Vertical  cross-section  and  plan  views  of  air  parcel  trajectories  near 
the  top  and  the  bottom  of  a  radially  symmetric  circulation. 

Fig.  8b.  Plan  view  of  air  parcel  trajectories  near  the  top  and  the  bottom  of  the 
same  radially  symmetric  circulation  which  is  combined  with  an  ambient  flow  and  a 
three-dimensional  model  of  fallout  particle  trajectories,  starting  at  points  P  and  R 
on  the  upper  boundary  of  the  circulation. 


deformed,  but  without  changes  of  magnitude  (continuity  principles  require  that 
a  stretching  in  the  y-direction  is  matched  by  shrinking  in  the  other  directions). 

In  such  cases,  the  expected  pattern  shape  and  location  should  be  estimated  from 
trajectory  computations. 

1.4  The  Influence  on  Fallout  of  Time-dependent  Vertical  Currents 

Both  the  time  of  onset  and  the  duration  of  fallout  at  a  point  on  the  ground 
can  be  influenced  by  time-dependent  vertical  circulations.  It  can  be  shown  that 
the  duration  of  fallout  at  a  point  on  the  ground  is  increased  beneath  intensifying 
updrafts  and  weakening  downdrafts,  and  decreased  beneath  weakening  updrafts 
and  intensifying  downdrafts.  Consider  particles  falling  through  an  updraft  which 
intensifies  with  time.  The  first  particle  falls  through  the  vertical  current  most 
rapidly,  and  the  last  most  slowly.  Therefore,  the  time  interval  between  the 
emergence  of  the  first  and  the  last  particle  is  greater  than  the  interval  separating 
their  entrance.  Consider  an  intensifying  downdraft.  The  first  particle  may  enter 
at  T^  and  emerge  at  T  +  15  minutes.  The  last  particle  may  enter  at  T^  and 
assisted  by  the  stronger  downdraft,  emerge  at  T2  +  10  minutes.  The  time  between 
emergences  at  the  bottom  is  5  minutes  less  than  the  interval  between  entrances  at 
the  top.  The  changed  duration  of  fallout  at  the  ground  is  reflected  linearly  in  a 
change  of  the  accumulation  of  fallout  at  the  ground. 

These  considerations  are  further  illustrated  by  Fig.  9  which  shows  schematic 
trajectories  of  particles  falling  through  an  intensifying  circulation.  In  the  rising 
currents ,  a  given  particle  experiences  stronger  low  level  horizontal  convergence 
for  a  longer  time  than  it  experiences  weaker  high  level  horizontal  divergence. 
Therefore,  the  points  of  exit  of  particles  converge  at  the  base  of  updraft  columns. 

A  similar  argument,  and  implications  of  continuity,  show  that  exit  points  must 
diverge  beneath  intensifying  downdrafts.  The  changed  spacing  of  trajectories 
at  the  base  of  the  circulation  is  associated  with  a  changed  vertical  distance 
separating  emerging  particles. 

These  considerations  are  still  further  illustrated  by  Fig.  10.  An  intensify¬ 
ing  circulation  and  fields  of  convergence  and  divergence  produce  a  deformation 
of  the  fallout  pattern.  Consistent  with  Fig.  9,  this  deformation  is  associated  with 
the  advection  of  sediment  from  the  downdraft  to  the  updraft  regions.  In  accordance 
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Fig  9  Schematic  indication  of  trajectories  (vertical)  and 
isochrones  (horizontal)  of  fallout  which  descends  through  an 
accelerating  local  circulation  (dashed  line). 
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Fig.  10.  Deformation  of  a  two-dimensional  atomic  cloud  by  an  intensifying 
local  circulation. 
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with  Eq.  (12),  this  advection  is  not  associated  with  changes  of  concentration  in 
any  particle  packets,  but  rather  with  stretching  and  shrinking  of  the  vertical 
extent  of  the  region  occupied  by  particles,  and  a  corresponding  change  in  the 
duration  of  fallout  at  the  ground.  Figure  11  shows  the  distribution  of  onset 
times  and  durations  of  fallout  at  the  ground  in  this  case. 

Interesting  limiting  cases  occur  when  the  intensifying  circulation  attains 
such  strength  that  the  maximum  updrafts  become  as  great  as  the  fall  speed  of 
particles.  Then  the  duration  of  fallout  beneath  the  downdraft  becomes  quite 
short,  while  beneath  the  updraft,  the  fallout,  once  started,  lasts  for  an  indefinitely 
long  time  over  an  indefinitely  small  horizontal  area  at  the  ground.  If  maximum 
updrafts  exceed  fall  speeds  before  some  particles  have  descended  past  the  point 
where  they  can  be  suspended,  the  onset  of  fallout  at  the  ground  at  the  point 
directly  beneath  the  updraft  core  will  be  indefinitely  delayed.  Of  course,  these 
examples  serve  a  pedagogic  purpose  principally,  since  circulation  changes  and 
diffusion  processes  will  always  alter  this  picture  greatly  in  the  real  world. 


1.5  Compressibility  of  the  Atmosphere  in  Relation  to  Fallout 

The  decrease  of  particle  fall  speed  and  the  increase  of  air  density  with 
decreasing  height  are  both  shown  by  Eq.  (11)  to  affect  the  concentration  of  fall¬ 
out.  Both  terms  on  the  right  of  Eq.  (11)  are  attributable  to  the  compressibility 
of  the  atmosphere. 

The  total  accumulation  F^  of  fallout  on  a  unit  area  on  the  ground  may  be 
defined  by 


T 

F  =  -  /  g  C  V  dt  (26) 

g  o  8  g 


where  C  is  the  mass  of  fallout  per  unit  volume  of  air  adjacent  to  the  ground,  V 
g  g 

is  the  velocity  of  descent  of  that  fallout,  and  T  is  the  length  of  time  that  fallout 

g 

continues  at  the  ground. 

It  is  readily  seen  that  T  is  unaffected  by  the  presence  of  a  local  steady  cir- 

g 

culaticn.  Consider  a  region  of  any  depth  (H  in  Fig.  12)  through  which  descends  a 
particle  packet  of  initial  depth  as  shown  in  Fig.  12.  If  successive  particles 
through  H  have  identical  velocity  histories,  the  time  for  the  trailing  edge  of  the 
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Fig.  11.  The  onset  time  (heavy  line)  and  duration  (thin  line)  of  fallout  at  the 
surface  below  the  intensifying  circulation  depicted  in  Fig.  10. 
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Fig.  12.  Vertical  shrinkage  of  a  two-dimensional  region  of  fallout  due  to  a 
decreasing  fall  velocity  of  the  particles. 
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column  of  depth  to  arrive  at  any  point  is  obviously  later  than  the  time  of 

arrival  of  the  leading  edge,  by  the  interval  -H^/V^.  When  the  packet  has 

descended  to  a  region  where  the  speed  of  descent  is  V  ,  its  height  H  is 

evidently  reduced  to  -  H^(V^/V^),  for  this  column,  descending  at  speed 

V_,  passes  a  point  in  the  same  time  as  the  column  of  depth  H  ,  descending  at 

1 

V 

The  decrease  of  fall  speed  following  descent  and  associated  shortening 
of  fallout  columns,  contributes  to  changed  concentrations  as  defined  by  Eq.  (11). 
This  process  alone,  however,  does  not  change  the  accumulation  of  fallout  at  the 
ground,  since,  in  Eq.  (26),  the  increase  of  C  would  be  exactly  compensated  by 
decreased  V  . 

g  j, 

When  a  vertical  current  is  present  in  the  column  of  depth  H,  other  con¬ 
centration-changing  factors  can  be  identified.  Since  the  fall  speed  decreases 
with  descent,  the  particle  packet  may,  for  example,  be  longer  exposed  to  hori¬ 
zontal  convergence  at  the  base  of  an  updraft  than  to  horizontal  divergence  at 
die  top.  The  net  effect  of  the  action  of  horizontal  convergence  on  descending 
particle  packets,  however,  is  reduced  in  a  compressible  atmosphere,  since  the 
level  of  no-horizontal  divergence  is  below  the  height  of  maximum  vertical  velocity. 
In  any  event,  the  change  of  horizontal  area  of  the  packet  is,  by  definition. 


1  dA  M  7*  3u  dv 
—  —  =  div  V  —  —  +  — . 
A  dt  2  dx  ay 


(27) 


The  fractional  change  of  cross-sectional  area  is  inversely  proportional  to  total 
deposit  of  fallout  per  unit  area,  i.e.; 

if  '  if  -  eXP  [  V“V2  7  dt]  ‘  eXP  ['  /«  diV2  ?  V^]  (28) 

Since  dt  =  ~dz/^r. +w)  (positive  for  downward  motions),  Eq.  (28)  can  be  used  to 
calculate  the  change  of  fallcwt  amount  when  V  ,  w,  and  3 In  p/dz  are  specified, 
with  div2  ^  defined  by  Eq.  (8a). 

Alternatively,  Eq.  (11)  can  be  used  to  calculate  the  concentration  at  the 
base  of  a  local  circulation;  this  multiplied  by  the  fall  speed  there,  V  ,  is  pro- 


portional  to  F  . 
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We  have 


C 

T  9V. 

T 

Inpr  -  - 
C2 

I  ~ 
0  92 

dt  +  / 

0 

C1 

*0  9V. 
C  t 

dz 

— -  =  exp 
C2 

H  dz 

V.+w 

l 

9  In  p 


dz  I 

Ywj 


Figures  13  and  14,  taken  from  [3],  show  variations  of  fall  speed  with  height 
for  various  diameters  of  fallout  particles ,  and  the  contribution  to  the  total  mass  of 
fallout  provided  by  particles  of  different  sizes.  Most  of  the  mass  is  typically 
contained  in  particles  about  60  p,  for  which 

.  -1  .  -  „ 


V,  =  V  (1  -  az)  m  sec 
i  0  ' 


-1  for  60p  particles) 


(a  *  0.25  X  10  ).  (31) 

For  an  example,  we  take  w  as  defined  by  Eq.  (13)  and,  in  the  troposphere, 

9  In  p  -4-1 

-  -  10  m  .  (32) 

Substitution  of  Eqs.  (31)  and  (12)  in  Eq.  (30)  gives  an  expression  which  can  be 
integrated  analytically,  but  this  is  extremely  cumbersome.  Much  time  and  labor 
can  be  saved  by  the  finite-difference  solution  given  below.  Since  {see  Eqs.  (27) 
and  (8)] 

1  dA  9u  9v  9w  9  In  p 


1  dA  3u  9v  _  9w  9  In 

A  dt  9x  9y  9z  9z 

we  find,  after  appropriate  substitutions,  that 

. . .-4 


I  dA 
A  dz 


(2z/H)  -  1 


v  ma-zM  t  L  _  s!) 

0  4w  1  Hi 

max  *  * 


Solutions  of  this  equation  are  plowed  in  Fig.  15. 

Figure  15  can  be  used  to  estimate  the  effect  of  vertical  circulations  on  fall¬ 
out  accumulations  attributable  to  the  compressibility  of  the  atmosphere.  It  should 
be  remembered  that  infinite  concentrations  are  always  confined  to  an  indefinitely 
small  area,  that  the  effect  of  an  updraft  to  increase  the  fallout  accumulation  is 
mitigated  by  its  corollary  action  in  delaying  the  arrival  of  fallout  at  the  surface, 
and  that  the  calculations  are  based  on  steady  circulations  without  diffusion.  Of 
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Fig.  13.  Terminal  velocities  of  spherical  particles 
as  a  function  of  altitude  and  particle  radius. 
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Fig.  14.  Size  distribution  of  radioactive  particles 
originating  from  the  mushroom  stem  of  an  atomic  cloud. 
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course,  the  total  amount  of  fallout  is  conserved  during  its  descent  through  the 
local  circulation. 

1.6  Local  Horizontal  Circulations  in  Relation  to  Fallout  Patterns 

Since  the  time  of  descent  of  particles  depends  on  the  sum  of  their  termi¬ 
nal  fall  speed  and  the  vertical  velocity  of  the  air,  the  pattern  of  fallout  which 

descends  through  a  nondiver  gent  local  circulation  +  =  0,  —  =  w  =  0) 

dx  ay  dz 

is  established  in  precisely  the  same  time  as  the  pattern  without  a  local  wind. 
However,  the  local  circulation  produces  a  horizontal  displacement  defined  by 

AX  =  TAu  =  -  |  Au,  (35) 

where  AX  is  the  displacement  due  to  the  local  circulation;  T  is  the  time  taken 
by  particles  to  descend  through  the  circulation  of  depth  Z,  Au  is  the  average 
difference,  at  a  descending  particle,  of  the  actual  air  velocity  and  the  velocity 
of  the  ambient  flow,  and  V  is  the  average  terminal  fall  speed.  When  the  average 
velocity  of  the  ambient  flow  is  relatively  small,  the  displacement  due  to  the  local 
circulation  dominates.  The  importance  of  the  local  circulation  is  proportional 
to  its  depth  and  intensity. 

To  illustrate  the  effect  of  a  local  horizontal  circulation  on  the  pattern  of 
fallout,  we  have  chosen  a  confluent -diffluent  wind  field  which  simulates  the  flow 
through  a  mountain  pass.  The  equations  of  motion  for  a  confluent-diffluent 
system  are 

u  =  uQ  ♦  ax;  v  =  vQ  -  ay  (36) 

in  which  a  is  positive  for  confluence  [10).  Each  air  trajectory  is  an  hyperbola 
associated  with  a  constant  k  such  that 

(uQ  +  az)  (vQ  -  ay)  =  k.  (37) 

-1  -3  -l 

In  this  study,  v^  =  0,  uQ  =  10  m  sec  ,  and  aQ  =  10  sec  .  The  resulting 
pattern  over  an  example  total  length  of  70  km  is  shown  in  Fig.  16. 

Figure  17a  shows  the  fallout  configuration  associated  with  detonation  of  a 
cylindrical  atomic  cloud  4  km  deep  and  5  km  in  diameter,  with  its  base  1  km 
above  the  ground.  Descent  of  particles  is  assumed  to  occur  wholly  in  the  environ¬ 
mental  wind  field.  Note  that  the  lightest  particles  are  most  widespread.  Figure 
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•  Confluent-di Shat  wind  field  sinnlatag  flow  through  a  mountain  pass.  The  b,  c,  and  d  give  the  posi- 

f  the  source  cloud  corresponding  with  the  fallout  patterns  B,  C,  and  D  in  Fig.  17. 
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^*8*  17.  Patterns  of  fallout  descending  from  a n  atomic  cloud  5  km  deep  through  a  wi zd  field  ooe  km  deep  (i 
'*  cwipifoi  with  that  in  undisturbed  parallel  flow:  A.  Lcdirtarbed  flow;  B.  Source  cloud  on  axis  of  sjmi 
€Sktnnre  of  the  circa  labor,  C.  Source  cloud  off  center;  and  D.  Scarce  eland  at  edge  of  circulation  near  fist 
km  high.  Part  of  the  fallout  lands  on  the  obstacle. 


17b  shows  the  pattern  due  to  an  on-axis  deto  nation  in  the  wind  field  of  Fig.  16; 
detonation  off  the  axis  at  point  c  in  Fig.  16  is  shown  in  Fig.  17c.  Pattern  d  shows 
another  off-center  explosion  close  to  the  boundary  of  an  obstacle  (mountain) 
with  a  flat  top  at  H  =  1  km.  A  part  of  the  material  is  translated  over  the  top 
with  an  undisturbed  velocity  u  and  settles.  The  remainder  is  caught  by  the  local 
circulation  through  which  it  descends  another  1000  m.  Note  that  the  influence  of 
the  local  circulation  in  all  cases  is  most  important  with  respect  to  the  smallest 
particles. 
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2.0  CONCLUSIONS 

This  study  has  considered  some  effects  of  local  circulation  on  fallout 
distributions.  The  role  of  precipitation  processes  and  diffusion  on  fallout  have 
been  omitted  from  consideration.  Principal  conclusions  follow. 

1.  When  horizontal  steady  winds  within  and  without  a  stabilized  atomic 
cloud  are  defined,  any  point  on  the  cloud  edge  can  be  associated  with  a  mass- 
per -unit-area  of  fallout  having  a  fall  speed  V. ,  a  point  on  the  ground,  and  the 
times  of  onset  and  termination  of  fallout  on  the  ground. 

2.  In  an  incompressible  atmosphere,  the  concentration  following  the 
motion  of  particles  with  constant  fall  speed  V  remains  constant.  In  the  com¬ 
pressible  atmosphere,  the  local  fallout  concentrations  beneath  updrafts  increase 
as  the  updraft  velocity  intensifies.  Changes  of  accumulation  by  a  factor  of  two 
or  more,  however,  depend  on  the  long  continued  settling  of  particles  in  vertical 
drafts  whose  speeds  are  nearly  the  same  as  the  terminal  fall  speeds  of  the 
particles.  Thus  large  local  increases  attributable  to  compressibility  occur 
only  under  highly  idealized  conditions  whose  occurrence  in  the  real  atmosphere 
should  be  investigated. 

3.  When  a  steady  local  circulation  occurs  about  a  horizontal  axis,  i.e., 
when  one  of  the  horizontal  wind  components  is  always  zero,  the  pattern  of  fall¬ 
out  eventually  established  is  identical  to  that  produced  by  the  ambient  non- 
divergent  flow  acting  alone.  However,  horizontal  winds  accompanying  vertical 
currents  which  are  symmetric  about  a  vertical  axis  redistribute  fallout  particles 
and  may  change  the  pattern  of  their  accumulation  at  the  ground.  In  either  case, 
the  time  of  arrival  of  particles  is  delayed  by  updrafts  and  expedited  by  down- 
drafts. 

4.  In  time-dependent  vertical  circulations,  the  duration  and  total  accumu¬ 
lation  of  fallout  with  constant  fall  velocity  is  increased  at  a  point  on  the  ground 
beneath  intensifying  updrafts  and  decreased  beneath  weakening  updrafts;  beneath 
downdrafts  the  converse  holds.  Large  changes  of  fallout  concentrations  occur 
in  association  with  small  changes  of  the  vertical  air  currents  when  their  speed 
is  nearly  the  same  as  the  fall  speed  of  particles.  In  just  these  cases  neglected 
diffusion  must  be  most  important,  since  the  time  for  descent  of  particles  is  long. 


34 


5.  The  pattern  of  fallout  which  descends  through  a  nondivergent  local 
circulation  is  established  in  precisely  the  same  time  as  the  pattern  without 
a  local  wind.  The  displacement  of  the  pattern  and  its  change  of  shape  by  the 
local  circulation  is  proportional  to  the  circulation’s  depth  and  intensity. 

6.  The  influence  of  the  local  circulation  is,  in  all  cases,  most  important 
with  respect  to  the  smallest  (slow  falling)  particles. 
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